Exploring the Genomic Landscape of Bacillus paranthracis PUMB_17 as a Proficient Phosphatidylcholine-Specific Phospholipase C Producer.

Phospholipases find versatile applications across industries, including detergent production, food modification, pharmaceuticals (especially in drug delivery systems), and cell signaling research. In this study, we present a strain of Bacillus paranthracis for the first time, demonstrating significant potential in the production of phosphatidylcholine-specific phospholipase C (PC-PLC). The investigation thoroughly examines the B. paranthracis PUMB_17 strain, focusing on the activity of PC-PLC and its purification process. Notably, the PUMB_17 strain displays extracellular PC-PLC production with high specific activity during the late exponential growth phase. To unravel the genetic makeup of PUMB_17, we employed nanopore-based whole-genome sequencing and subsequently conducted a detailed genome annotation. The genome comprises a solitary circular chromosome spanning 5,250,970 bp, featuring a guanine-cytosine ratio of 35.49. Additionally, two plasmids of sizes 64,250 bp and 5845 bp were identified. The annotation analysis reveals the presence of 5328 genes, encompassing 5186 protein-coding sequences, and 142 RNA genes, including 39 rRNAs, 103 tRNAs, and 5 ncRNAs. The aim of this study was to make a comprehensive genomic exploration that promises to enhance our understanding of the previously understudied and recently documented capabilities of Bacillus paranthracis and to shed light on a potential use of the strain in the industrial production of PC-PLC.

Illuminating the Genomic Landscape of Lactiplantibacillus plantarum PU3-A Novel Probiotic Strain Isolated from Human Breast Milk, Explored through Nanopore Sequencing.

Lactiplantibacillus plantarum stands out as a remarkably diverse species of lactic acid bacteria, occupying a myriad of ecological niches. Particularly noteworthy is its presence in human breast milk, which can serve as a reservoir of probiotic bacteria, contributing significantly to the establishment and constitution of infant gut microbiota. In light of this, our study attempted to conduct an initial investigation encompassing both genomic and phenotypic aspects of the L. plantarum PU3 strain, that holds potential as a probiotic agent. By employing the cutting-edge third-generation Nanopore sequencing technology, L. plantarum PU3 revealed a circular chromosome of 3,180,940 bp and nine plasmids of various lengths. The L. plantarum PU3 genome has a total of 2962 protein-coding and non-coding genes. Our in-depth investigations revealed more than 150 probiotic gene markers that unfold the genetic determinants for acid tolerance, bile resistance, adhesion, and oxidative and osmotic stress. The in vivo analysis showed the strain's proficiency in utilizing various carbohydrates as growth substrates, complementing the in silico analysis of the genes involved in metabolic pathways. Notably, the strain demonstrated a pronounced affinity for D-sorbitol, D-mannitol, and D-Gluconic acid, among other carbohydrate sources. The in vitro experimental verification of acid, osmotic and bile tolerance validated the robustness of the strain in challenging environments. Encouragingly, no virulence factors were detected in the genome of PU3, suggesting its safety profile. In search of beneficial properties, we found potential bacteriocin biosynthesis clusters, suggesting its capability for antimicrobial activity. The characteristics exhibited by L. plantarum PU3 pave the way for promising strain potential, warranting further investigations to unlock its full capacity and contributions to probiotic and therapeutic avenues.

The Fascinating World of Plant Non-Coding RNAs.

Non-coding RNAs (ncRNAs), such as microRNAs (miRNAs), short interfering RNAs (siRNAs), and long non-coding RNAs (lncRNAs), have emerged as pivotal regulators within the plant kingdom [...].

Liquid Biopsy: Current Status and Future Perspectives.

Since the discovery of the Bence Jones protein in the middle to late 1800s and the subsequent identification of the carcinoembryonic antigen and alpha-fetoprotein in the 1970s, it has been demonstrated that the analysis of biofluids is essential to the diagnostic and follow-up processes of cancer [...].

Identification of Bulgarian Sourdough Microbiota by Metagenomic Approach Using Three Commercially Available DNA Extraction Protocols.

Research background: Sourdough is a spontaneously formed, complex microbial ecosystem of various lactic acid bacteria (LAB) and yeast which, by producing specific metabolites, determines the quality of the baked products. In order to design and control the sourdough with preferred nutritional characteristics, it is crucial that the LAB diversity of the product of interest be elucidated. Experimental approach: Using the opportunities of next-generation sequencing (NGS) of the V1-V3 hypervariable gene region of 16S rRNA, we studied the microbial ecosystem of a whole grain sourdough made of Triticum monococcum, originating from Southwestern Bulgaria. Since the DNA extraction method is considered crucial for the accuracy of the sequencing results, as it can introduce significant differences in the examined microbiota, we used three different commercial kits for DNA isolation and analyzed their impact on the observed bacterial diversity. Results and conclusions: All three DNA extraction kits provided bacterial DNA which passed quality control and was successfully sequenced on Illumina MiSeq platform. The results received from the different DNA protocols showed variations in the microbial profiles. Alpha diversity indices (ACE, Chao1, Shannon, and Simpson) were also different among the three groups of results. Nevertheless, a strong dominance of phylum Firmicutes, class Bacilli, order Lactobacillales, represented mostly by family Lactobacillaceae, genus Lactobacillus (relative abundance of 63.11-82.28%) and family Leuconostocaceae, genus Weissella (relative abundance of 3.67-36.31%) was observed. Lactiplantibacillus plantarum and Levilactobacillus brevis with relative abundance of 16.15-31.24% and 6.21-16.29% respectively, were the two dominant species identified in all three DNA isolates. Novelty and scientific contribution: The presented results give insight into the taxonomic composition of bacterial community of a specific Bulgarian sourdough. Having in mind that the sourdough is a difficult matrix for DNA isolation on the one hand, and that there is no standardized DNA extraction protocol for this matrix on the other hand, this pilot study aims to give a small contribution to the future establishment and validation of such a protocol, which will allow accurate assessment of the specific microbiota of sourdough samples.

16S-rRNA-Based Metagenomic Profiling of the Bacterial Communities in Traditional Bulgarian Sourdoughs.

Sourdoughs (SDs) are spontaneously formed microbial ecosystems composed of various species of lactic acid bacteria (LAB) and acid-tolerant yeasts in food matrices of cereal flours mixed with water. To date, more than 90 LAB species have been isolated, significantly impacting the organoleptic characteristics, shelf life, and health properties of bakery products. To learn more about the unique bacterial communities involved in creating regional Bulgarian sourdoughs, we examined the metacommunities of five sourdoughs produced by spontaneous fermentation and maintained by backslopping in bakeries from three geographic locations. The 16S rRNA gene amplicon sequencing showed that the former genus Lactobacillus was predominant in the studied sourdoughs (51.0-78.9%). Weissella (0.9-42.8%), Herbaspirillum (1.6-3.8%), Serratia (0.1-11.7%), Pediococcus (0.2-7.5%), Bacteroides (0.1-1.3%), and Sphingomonas (0.1-0.5%) were also found in all 5 samples. Genera Leuconostoc, Enterococcus, Bacillus, and Asaia were sample-specific. It is interesting to note that the genus Weissella was more abundant in wholegrain samples. The greatest diversity at the species level was found in the former genus Lactobacillus, presented in the sourdough samples with 13 species. The UPGMA cluster analysis clearly demonstrated similarity in species' relative abundance between samples from the same location. In addition, we can conclude that the presence of two main clusters-one including samples from mountainous places (the cities of Smolyan and Bansko) and the other including samples from the city of Ruse (the banks of the Danube River)-may indicate the impact of climate and geographic location (e.g., terrain, elevation, land use, and nearby water bodies and their streams) on the abundance of microbiome taxa. As the bacterial population is crucial for bread standardization, we expect the local bakery sector to be interested in the relationship between process variables and their effect on bacterial dynamics described in this research study.

The Multiverse of Plant Small RNAs: How Can We Explore It?

Plant small RNAs (sRNAs) are a heterogeneous group of noncoding RNAs with a length of 20-24 nucleotides that are widely studied due to their importance as major regulators in various biological processes. sRNAs are divided into two main classes-microRNAs (miRNAs) and small interfering RNAs (siRNAs)-which differ in their biogenesis and functional pathways. Their identification and enrichment with new structural variants would not be possible without the use of various high-throughput sequencing (NGS) techniques, allowing for the detection of the total population of sRNAs in plants. Classifying sRNAs and predicting their functional role based on such high-performance datasets is a nontrivial bioinformatics task, as plants can generate millions of sRNAs from a variety of biosynthetic pathways. Over the years, many computing tools have been developed to meet this challenge. Here, we review more than 35 tools developed specifically for plant sRNAs over the past few years and explore some of their basic algorithms for performing tasks related to predicting, identifying, categorizing, and quantifying individual sRNAs in plant samples, as well as visualizing the results of these analyzes. We believe that this review will be practical for biologists who want to analyze their plant sRNA datasets but are overwhelmed by the number of tools available, thus answering the basic question of how to choose the right one for a particular study.

The Dynamics of miR-449a/c Expression during Uterine Cycles Are Associated with Endometrial Development.

The human endometrium is a highly dynamic tissue. Increasing evidence has shown that microRNAs (miRs) play essential roles in human endometrium development. Our previous assay, based on small RNA-sequencing (sRNA-seq) indicated the complexity and dynamics of numerous sequence variants of miRs (isomiRs) that can act together to control genes of functional relevance to the receptive endometrium (RE). Here, we used a greater average depth of sRNA-seq to detect poorly expressed small RNAs. The sequencing data confirmed the up-regulation of miR-449c and uncovered other members of the miR-449 family up-regulated in RE-among them miR-449a, as well as several isoforms of both miR-449a and miR-449c, while the third family member, miR-449b, was not identified. Stem-looped RT-qPCR analysis of miR expression at four-time points of the endometrial cycle verified the increased expression of the miR-449a/c family members in RE, among which the 5' isoform of miR-449c-miR-449c.1 was the most strongly up-regulated. Moreover, we found in a case study that the expression of miR-449c.1 and its precursor correlated with the histological assessment of the endometrial phase and patient age. We believe this study will promote the clinical investigation and application of the miR-449 family in the diagnosis and prognosis of human reproductive diseases.

Transcriptome Analysis Reveals Dynamic Cultivar-Dependent Patterns of Gene Expression in Potato Spindle Tuber Viroid-Infected Pepper.

Potato spindle tuber viroid (PSTVd) infects various plants. PSTVd pathogenesis is associated with interference with the cellular metabolism and defense signaling pathways via direct interaction with host factors or via the transcriptional or post-transcriptional modulation of gene expression. To better understand host defense mechanisms to PSTVd infection, we analyzed the gene expression in two pepper cultivars, Capsicum annuum Kurtovska kapia (KK) and Djulunska shipka (DS), which exhibit mild symptoms of PSTVd infection. Deep sequencing-based transcriptome analysis revealed differential gene expression upon infection, with some genes displaying contrasting expression patterns in KK and DS plants. More genes were downregulated in DS plants upon infection than in KK plants, which could underlie the more severe symptoms seen in DS plants. Gene ontology enrichment analysis revealed that most of the downregulated differentially expressed genes in both cultivars were enriched in the gene ontology term photosynthesis. The genes upregulated in DS plants fell in the biological process of gene ontology term defense response. We validated the expression of six overlapping differentially expressed genes that are involved in photosynthesis, plant hormone signaling, and defense pathways by quantitative polymerase chain reaction. The observed differences in the responses of the two cultivars to PSTVd infection expand the understanding of the fine-tuning of plant gene expression that is needed to overcome the infection.

Coupling miR/isomiR and mRNA Expression Signatures Unveils New Molecular Layers of Endometrial Receptivity.

Embryo implantation depends on endometrial receptivity (ER). To achieve ER, the preparation of the uterine lining requires controlled priming by ovarian hormones and the expression of numerous genes in the endometrial tissue. microRNAs (miRs) have emerged as critical genetic regulators of ER in fertility and of the diseases that are associated with infertility. With the rapid development of next-generation sequencing technologies, it has become clear that miR genes can produce canonical miRs and variants-isomiRs. Here, we describe miR/isomiR expression dynamics across the four time points of natural chorionic gonadotropin (hCG)-administered cycles. Sequencing of the small RNAs (sRNA-seq) revealed that the most significant expression changes during the transition from the pre-receptive to the receptive phase occurred in the isomiR families of miR-125a, miR-125b, miR-10a, miR-10b, miR-449c, miR-92a, miR-92b, and miR-99a. Pairing the analysis of the differentially expressed (DE) miRs/isomiRs and their predicted DE mRNA targets uncovered 280 negatively correlating pairs. In the receptive endometrium, the 5'3'-isomiRs of miR-449c, which were among the most highly up-regulated isomiRs, showed a negative correlation with their target, transcription factor (TF) MYCN, which was down-regulated. Joint analysis of the miR/isomiR and TF expression identified several regulatory interactions. Based on these data, a regulatory TF-miR/isomiR gene-target circuit including let7g-5p and miR-345; the isomiR families of miR-10a, miR-10b, miR-92a, and miR-449c; and MYCN and TWIST1 was proposed to play a key role in the establishment of ER. Our work uncovers the complexity and dynamics of the endometrial isomiRs that can act cooperatively with miRs to control the functionally important genes that are critical to ER. Further studies of miR/isomiR expression patterns that are paired with those of their target mRNAs may provide a more in-depth picture of the endometrial pathologies that are associated with implantation failure.

miRGalaxy: Galaxy-Based Framework for Interactive Analysis of microRNA and isomiR Sequencing Data.

Tools for microRNA (miR) sequencing data analyses are broadly used in biomedical research. However, the complexity of computational approaches still remains a challenge for biologists with scarce experience in data analytics and bioinformatics. Here, we present miRGalaxy, a Galaxy-based framework for comprehensive analysis of miRs and their sequence variants-miR isoforms (isomiRs). Though isomiRs are commonly reported in deep-sequencing experiments, their detailed structure complexity and specific differential expression (DE) remain not fully examined by the majority of the available analysis tools. miRGalaxy encompasses biologist-user-friendly tools and workflows dedicated to the analysis of the isomiR-ome and its complex behavior in various biological samples. miRGalaxy is developed as a modular, accessible, redistributable, shareable, and user-friendly framework for scientists working with small RNA (sRNA)-seq data. Due to its modular workflow, advanced users can customize the steps and tools for their needs. In addition, the framework provides an analysis report where the significant output results are summarized in charts and visualizations. miRGalaxy can be accessed via preconfigured Docker image flavor and a Toolshed installation if the user already has a running Galaxy instance. Over the last decade, studies on the expression of miRs and isomiRs in normal and deregulated tissues have led to the discovery of their potential as diagnostic biomarkers. The detection of miRs in biofluids further expanded the exploration of the miR repertoire as a source of liquid biopsy biomarkers. Here we show the miRGalaxy framework application for in-depth analysis of the sRNA-seq data from two different biofluids, milk and plasma, to identify, annotate, and discover specific differentially expressed miRs and isomiRs.

Analysis of the complete mitochondrial genome sequence of the resurrection plant Haberlea rhodopensis.

Haberlea rhodopensis is a paleolithic tertiary relict species that belongs to the unique group of resurrection plants sharing remarkable tolerance to desiccation. When exposed to severe drought stress, this species shows an ability to maintain structural integrity of its deactivated photosynthetic apparatus, which easily reactivates upon rehydration. In addition to its homoiochlorophyllous nature, the resurrection capability of H. rhodopensis is of particular importance to the global climate change mitigation. In this study, we sequenced, assembled, and analyzed the mitochondrial (mt) genome of H. rhodopensis for the first time. The master circle has a typical circular structure of 484 138 bp in length with a 44.1% GC content in total. The mt genome of H. rhodopensis contains 59 genes in total, including 35 protein-coding, 21 tRNAs, and 3 rRNAs genes. 7 tandem repeats and 85 simple sequence repeats (SSRs) are distributed throughout the mt genome. The alignment of 20 plant mt genomes confirms the phylogenetic position of H. rhodopensis in the Lamiales order. Our comprehensive analysis of the complete mt genome of H. rhodopensis is a significant addition to the limited database of organelle genomes of resurrection species. Comparative and phylogenetic analysis provides valuable information for a better understanding of mitochondrial molecular evolution in plants.

isomiRs-Hidden Soldiers in the miRNA Regulatory Army, and How to Find Them?

Numerous studies on microRNAs (miRNA) in cancer and other diseases have been accompanied by diverse computational approaches and experimental methods to predict and validate miRNA biological and clinical significance as easily accessible disease biomarkers. In recent years, the application of the next-generation deep sequencing for the analysis and discovery of novel RNA biomarkers has clearly shown an expanding repertoire of diverse sequence variants of mature miRNAs, or isomiRs, resulting from alternative post-transcriptional processing events, and affected by (patho)physiological changes, population origin, individual's gender, and age. Here, we provide an in-depth overview of currently available bioinformatics approaches for the detection and visualization of both mature miRNA and cognate isomiR sequences. An attempt has been made to present in a systematic way the advantages and downsides of in silico approaches in terms of their sensitivity and accuracy performance, as well as used methods, workflows, and processing steps, and end output dataset overlapping issues. The focus is given to the challenges and pitfalls of isomiR expression analysis. Specifically, we address the availability of tools enabling research without extensive bioinformatics background to explore this fascinating corner of the small RNAome universe that may facilitate the discovery of new and more reliable disease biomarkers.

Transposon-associated polymorphisms of stress-responsive gene promoters in selected accessions of Arabidopsis thaliana.

Genetic diversity caused by transposable element movement can play an important role in plant adaptation to local environments. Regarding genes, transposon-induced alleles were mostly related to gene bodies and a few of them to promoter regions. In this study, promoter regions of 9 stress-related genes were searched for transposable element insertions in 12 natural accessions of Arabidopsis thaliana. The promoter screening was performed via PCR amplification with primers designed to flank transposable element insertions in the promoter regions of the reference accession Col-0. Transposable element-associated insertion/deletion (indel) polymorphisms were identified in 7 of the 12 promoter loci across studied accessions that can be developed further as molecular markers. The transposable element absence in the promoter regions of orthologous genes in A. lyrata indicated that the insertion of these transposable elements in A. thaliana lineage had occurred after its divergence from A. lyrata. Sequence analysis of the promoter regions of CML41 (Calmodulin-like protein 41) and CHAP (chaperone protein dnaJ-related) confirmed the indel polymorphic sites in four accessions - Col-0, Wassilewskija, Shahdara, and Pirin. The observed indel polymorphism of the CHAP promoter region was associated with specific gene expression profiles in the different accessions grown at a normal and elevated temperature in a plant growth chamber. The collected data can be a starting point for gene expression profiling studies under conditions resembling the natural habitats of accessions.

Metagenomic profiling of the microbial freshwater communities in two Bulgarian reservoirs.

Microorganisms inhabiting freshwater environments are an integral part of the aquatic ecosystems. Very few data are available regarding the profiles of the microbial communities in the reservoirs in Bulgaria, despite their key role in the biogeochemical processes. In the present study, we provide the first comprehensive metagenomic analysis on the planktonic bacterial diversity of two large and economically important Bulgarian reservoirs (Batak and Tsankov Kamak) using next-generation sequencing of 16S ribosomal RNA gene (16S rRNA). Analysis of the metagenomic amplicon datasets, including quality filtering, clustering of Operational Taxonomic Units and taxonomy assignment revealed that 78.45% of the microbial communities between the two reservoirs were overlapping. The diversity (H) and Pielou's evenness (J) indices declined along the longitudinal axis of both reservoirs. The estimated values for the Shannon diversity index are typically observed in oligotrophic lakes. The microbial communities of both reservoirs were dominated by Proteobacteria, followed by Actinobacteria and Bacteroidetes all comprised over 95% of the relative abundance, regardless of the reservoir's large hydrogeological differences. The bacterioplankton was characterized by high phylogenetic heterogeneity in the taxonomic structure, being distributed among 211 genera. The genera Limnohabitans and Rhodoferax held the absolute predominance, implying their significance in the aquatic food webs. The obtained data can contribute to the better systematic understanding of the microbial diversity of freshwater environments.

Chloroplast Genome Analysis of Resurrection Tertiary Relict Haberlea rhodopensis Highlights Genes Important for Desiccation Stress Response.

Haberlea rhodopensis is a paleolithic tertiary relict species, best known as a resurrection plant with remarkable tolerance to desiccation. When exposed to severe drought stress, H. rhodopensis shows an ability to maintain the structural integrity of its photosynthetic apparatus, which re-activates easily upon rehydration. We present here the results from the assembly and annotation of the chloroplast (cp) genome of H. rhodopensis, which was further subjected to comparative analysis with the cp genomes of closely related species. H. rhodopensis showed a cp genome size of 153,099 bp, harboring a pair of inverted repeats (IR) of 25,415 bp separated by small and large copy regions (SSC and LSC) of 17,826 and 84,443 bp. The genome structure, gene order, GC content and codon usage are similar to those of the typical angiosperm cp genomes. The genome hosts 137 genes representing 70.66% of the plastome, which includes 86 protein-coding genes, 36 tRNAs, and 4 rRNAs. A comparative plastome analysis with other closely related Lamiales members revealed conserved gene order in the IR and LSC/SSC regions. A phylogenetic analysis based on protein-coding genes from 33 species defines this species as belonging to the Gesneriaceae family. From an evolutionary point of view, a site-specific selection analysis detected positively selected sites in 17 genes, most of which are involved in photosynthesis (e.g., rbcL, ndhF, accD, atpE, etc.). The observed codon substitutions may be interpreted as being a consequence of molecular adaptation to drought stress, which ensures an evolutionary advantage to H. rhodopensis.

ChloroMitoSSRDB 2.00: more genomes, more repeats, unifying SSRs search patterns and on-the-fly repeat detection.

Organelle genomes evolve rapidly as compared with nuclear genomes and have been widely used for developing microsatellites or simple sequence repeats (SSRs) markers for delineating phylogenomics. In our previous reports, we have established the largest repository of organelle SSRs, ChloroMitoSSRDB, which provides access to 2161 organelle genomes (1982 mitochondrial and 179 chloroplast genomes) with a total of 5838 perfect chloroplast SSRs, 37 297 imperfect chloroplast SSRs, 5898 perfect mitochondrial SSRs and 50 355 imperfect mitochondrial SSRs across organelle genomes. In the present research, we have updated ChloroMitoSSRDB by systematically analyzing and adding additional 191 chloroplast and 2102 mitochondrial genomes. With the recent update, ChloroMitoSSRDB 2.00 provides access to a total of 4454 organelle genomes displaying a total of 40 653 IMEx Perfect SSRs (11 802 Chloroplast Perfect SSRs and 28 851 Mitochondria Perfect SSRs), 275 981 IMEx Imperfect SSRs (78 972 Chloroplast Imperfect SSRs and 197 009 Mitochondria Imperfect SSRs), 35 250 MISA (MIcroSAtellite identification tool) Perfect SSRs and 3211 MISA Compound SSRs and associated information such as location of the repeats (coding and non-coding), size of repeat, motif and length polymorphism, and primer pairs. Additionally, we have integrated and made available several in silico SSRs mining tools through a unified web-portal for in silico repeat mining for assembled organelle genomes and from next generation sequencing reads. ChloroMitoSSRDB 2.00 allows the end user to perform multiple SSRs searches and easy browsing through the SSRs using two repeat algorithms and provide primer pair information for identified SSRs for evolutionary genomics.Database URL: http://www.mcr.org.in/chloromitossrdb.

High-temperature effect on genes engaged in DNA methylation and affected by DNA methylation in Arabidopsis.

Along with its essential role in the maintenance of genome integrity, DNA methylation takes part in regulation of genes which are important for plant development and stress response. In plants, DNA methylation process can be directed by small RNAs in process known as RNA-directed DNA methylation (RdDM) involving two plant-specific RNA polymerases - PolIV and PolV. The aim of the present study was to investigate the effect of heat stress on the expression of genes encoding key players in DNA methylation - DNA methyltransferase (MET1, CMT3, and DRM2), the largest subunits of PoIIV and PolV (NRPD1 and NRPE1 respectively) and the DNA demethylase ROS1. We also examined the high-temperature effect on two protein-coding genes - At3g50770 and At5g43260 whose promoters contain transposon insertions and are affected by DNA-methylation, as well as on the AtSN1, a SINE-like retrotransposon. To assess the involvement of PolIV and PolV in heat stress response, the promoter methylation status and transcript levels of these genes were compared between wild type and double mutant lacking NRPD1 and NRPE1. The results demonstrate coordinated up-regulation of the DRM2, NRPD1 and NRPE1 in response to high temperature and suggest that PolIV and/or PolV might be required for the induction of DRM2 expression under heat stress. The ROS1 expression was confirmed to be suppressed in the mutant lacking active PolIV and PolV that might be a consequence of abolished DNA methylation. The increased expression of At3g50770 in response to elevated temperature correlated with reduced promoter DNA methylation, while the stress response of At5g43260 did not show inverse correlation between promoter methylation and gene expression. Our results also imply that PolIV and/or PolV could regulate gene expression under stress conditions not only through RdDM but also by acting in other regulatory processes.

Insight into small RNA abundance and expression in high- and low-temperature stress response using deep sequencing in Arabidopsis.

Small RNA profiling and assessing its dependence on changing environmental factors have expanded our understanding of the transcriptional and post-transcriptional regulation of plant stress responses. Insufficient data have been documented earlier to depict the profiling of small RNA classes in temperature-associated stress which has a wide implication for climate change biology. In the present study, we report a comparative assessment of the genome-wide profiling of small RNAs in Arabidopsis thaliana using two conditional responses, induced by high- and low-temperature. Genome-wide profiling of small RNAs revealed an abundance of 21 nt small RNAs at low temperature, while high temperature showed an abundance of 21 nt and 24 nt small RNAs. The two temperature treatments altered the expression of a specific subset of mature miRNAs and displayed differential expression of a number of miRNA isoforms (isomiRs). Comparative analysis demonstrated that a large number of protein-coding genes can give rise to differentially expressed small RNAs following temperature shifts. Low temperature caused accumulation of small RNAs, corresponding to the sense strand of a number of cold-responsive genes. In contrast, high temperature stimulated the production of small RNAs of both polarities from genes encoding functionally diverse proteins.